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I. IN TROD OCTIOH 

Guided missiles are classified into four broad catego- 
ries, depending on launch and target position characteris- 
tics. These categories are (1) . air -to - air (2). air - 
to - ground (3) . surface - to - air and (4) . surface -to - 
surface. Each category of the above will employ one or more 
of the following guidance schemas; programmed command, 
liae-of-sight, lead-angle, proportional navigation homing 
and inertial. The beam rider guidance is included in the 
line-of-sight guidance. A number of missiles also use a 
combination of these methods. For example, the initial part 
of the missile trajectory may use programmed guidance while 
the terminal phase may use beam-rider. 

This thesis discusses the surface-to-air missile 
controlled by on-off, thrust vector, control. Consideration 
was given to determine the effects of the two-level relay 
and the saturation linear control. In order to verify the 
results, it was tested by using the type of control for 
three different types of missile-target scenarios: 

(1) . LOS command against non-maneuvering target 

(2) . LOS command against maneuvering target 

(3) . Pseudo-LOS command against non-maneuvering target. 

In chapter 2, a discussion of a line-of-sight guidance 

was presented and a practical example of it was shown in 
chapter 3. The general concept of on-off control was 
described in chapter 4. The simulation results of the basic 
command to line-of-sight against both a non- maneuvering and 
a maneuvering target were shown in chapters 5 and that of 
pseudo- LOS case was in chapter 6. Finally, a discussion of 
two-level relay and saturating linear control was presented 
in chapter 7. A table of variables which were used in this 



9 



was 



thesis is shewn 
accomplished us 



in the Appendix A. Computer simul; 
ng Digital Simulation Language, DSL 



•cion 



10 



II. OVEHVIEW OF LINE::OF:;SIGHr GUIDANCE CONTROL 

A LOS system can be called a "3-point” guidance system 
since there is one point which defines the tracker, another 
the target and a third which defines the position of the 
missile. The object of the guidance system is to constrain 
the missile to lie as nearly as possible on the line joining 
the tracker and the target called the Line Of Sight (LOS). 
The concept is simple and can be implemented in many ways; 
perhaps it is this apparent simplicity which explains why 
many of the guided weapon systems as yet designed are LOS 
system. 

Consider a target flying straight and at constant speed, 
and a missile flying at a different angle but constant 
speed, having been launched when the target occupies a posi- 
tion TO (see Figure 2.1). 

After intervals of time of 1,2,3 etc seconds the LOS is 
shown as OT1,OT2,OT3 etc. Since the missile ideally always 
lies on these lines the flight path will be a curved one, 
for an approaching target, the curvature becomes increas- 
ingly severe towards the end of the engagement. We note that 
the tangent to the flight path at any one point defines the 
instantaneous direction of the missile velocity. It is seen 
that the missile velocity vector will, in general, not be 
directed along the LOS; towards the end of the engagement it 
may be at a considerable angle to it [Ref. 1]. 

In an actual situation the guidance signals transmitted 
to the missile are the demanded lateral accelerations 
(LATAX) in two axes at the right angles to the beam. These 
demands are resolved into missile axes within the missile. 
An error compensation term endeavouring to keep the error 
off the beam ( equal to zero. 

1 1 




Figare 2.1 aissile Target Encounter with LOS Guidance 

A basic geometry and a simplified guidance loop are 
shown in Figure 2.2 and Figure 2.3. 

Suppose that the cross range error (CRE) of Figure 2.2 
can be measured either directly or by means of the angular 
difference between OT and oa, together with some knowledge 
of missile range (Rm) , then 

CRE = -<5'm ) (2. D 

If this error off the beam is used as an acceleration demand 
0, it needs some damping so that good response characteris- 
tics are obtained. A dynamic equation of the form 

CRE = G1 (CRE) + G2 (CRE) (2.2) 

needs to be satisfied, where G1 and G2 are constants. This 
necessity leads immediately to the consideration of a 
filtered error. In the presence of noise on the sight-line. 
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Figure 2,2 Basic Seometry 




Figure 2.3 Simplified Guidance Loop of LOS Guidance 
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and hence on the cross range error, CRS, such a filrar 
design is not simple and becomes a compromise between 
requirements for smoothing the noisa and giving an adequate 
response to a demand. Modern techniques allow filters to be 
designed statistically if some knowledge of the noise char- 
acteristics is available or can be assumed. Figure 2.3 shows 
the position of such a filter F(s) in the guidance loop. It 
includes a gain G, and the acceleration demand is 

U = F (s) Rm ( (5t - ) (2.3) 

The missile transfer function is represented by A(s) and 
when the achieved acceleration is doubly integrated and 
divided by Rm it represents a new measure of rhe missile 
beam angle (d ^) , thus closing the loop when differenced with 
the target beam angle (dt) • 

While this concept is simply a LOS or beam riding 
guidance situation it is by no means as clear in homing how 
a guidance law can be devised in the absence of information 
on missile and target positions [Ref. 2]. 
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III. TYPI CAL engagement OQ2M£e 



In order to provide a "vehicle" through which to better 
understand the basic aspects of command to line-of-sight 
guidance methodology, the engagement sequence of a short- 
range, air-defense, missile system is described. The Soland 
system was selected because the general operational aspects 
of the system are available at the unclassified level 
[Bef. 4]. 

The entry of one or more aerial targets into the range 
of the search radar is indicated to the Roland vehicle 
commander by an audible tone. At the sane time, a synthetic 
display of the targets appears on a screen to give the 
commmander all the information needed to select the most 
threatening target. The screen images are different for 
friendly and enemy targets. Also, the entry of the target 
into the missile envelope, utilizing target advanced- range 
computations, is indicated by a change in the display. With 
the search antenna raised and the search radar activated, 
target acquisition is possible even when the vehicle is in 
motion. 

There are three modes of identification, friend or foe 
(IFF) interrogation: automatic, manual, and automatic 

within a given range. 

When the commander has recognized a target as hostile 
and decided to engage it, he places a cursor over the screen 
image.. This automatically brings the turret to bear and 
tracking can commence in either the "radar" or "optical" 
modes. 

In the "radar" mode, the tracking radar automatically 
accepts target designation from the search radar, searches 
for, locks onto, and tracks the target. 
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In the ''optical" mode, the aimer searches for the target 
in elevation with an optical sight. To aid him an elecrronic 
instrument displays the maximum theoretical elevation for 
the search. When the aimer has acquired the target in his 
cross-hairs, he keeps the target in his sight by manipu- 
lating a control stick. This control keeps zhe target prop- 
erly positioned by moving the turret in azimuth and 
swivelling a mirror in elevation. 

As soon as the commander confirms that the target is 
within missile range, he initiates the firing sequence in 
the "radar” mode, or authorizes "optical" mode firing 

through a command displayed in the aimer's sight. The aimer, 
then, can initiate the firing sequence. 

The missile is guided by a command to line-of-sight 
technique. This means that the target is tracked optically 
or by radar and the deviation of the missile from this line 
of sight is determined and corrected by a guidance command. 
The commander may switch from "radar" to "optical" and back 
again, as desired, even after the missile has been launched. 

Target tracking and determination of the missile's devi- 
ation from the line of sight are different for each mode. In 
the "radar" mode, the guidance radar has two receiving 
channels. One is used for target tracking and the other is 
used to locate the missile in xhe radar lobe through recep- 
tion of the missile's radio frequency beacons. By comparing 
these angles, an error between the missile and the target 
line of sight can be determined. In the "oprical" mode, a 
biaxially-stabilized mirror is manually controlled to keep 
the target vertically in the aimer's sight and the turret is 
rotated to the azimuth of the target line of sight. An 
infrared goniometer is mounted to provide misile angle from 
the tracker by following flares mounted on the rear of the 
missile. Then, a deviation of the missile angle from the 
target line of sight can be determined. 
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Twc groups of signals are introduced into the command 
computer; the velocity of the line of sight in azimuth and 
elevation, and the deviation of the missile from the line of 
sight in azimuth and elevation. Based upon data from the 
line-of-sight movement and the angular deviations of the 
missile, the necessary guidance signals are calculated. 

The guidance signals are relayed to the missile by a 
transmitter with highly directional characteristics. The 
command-transmitting antenna is slaved to the missile angle 
in both azimuth and elevation. It, therefore, is trained on 
the missile continuously. 

The side forces required for missile course corrections 
are produced through deflection of the exhaust jet of the 
sustainer motor by spoilers at the rear of the missile 
(thrust- vector control). 

When the missile reaches the point of impact with the 
target, the warhead is detonated by either percussion, 
contact fuse or the radio- frequency , proximity fuse. The 
warhead consists of 'a radial-effect, multiple-fragmentation 
charge. 

Figure 3,1 presents an operational schematic of the 
basic Roland missile system operation. 

The computer simulations contained herein are generic in 
nature within the command to line-of-sight guided-missile 
type and have only reasonable estimates of missile capabil- 
ities introduced. This ensures unclassified results. At the 
same time, the simulations are of sufficient complexity to 
properly weigh the relative merits of the guidance 
variations discussed [Ref, 3]. 
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Figure 3. 1 



Boland Hissile Systea 



Operational Schematic 
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IV. O M-OFF (BAN^BANG) CONTROL 



As discussed before, LOS guidance maintains a missile 
position on the LOS, Usually missile position has a cross 
range error (CRE) and we want to reduce this error to zero 
in the minimum time. This kind of problem can be solved by 
using the on-off control. The basic concept of this is that; 



Given a system for which the drive is limited (has a 
maximum or saturation value) , the fastssr response is 
obtained if maximum forward drive is applied at t = 0, 
and is reversed at 'a proper instant t = t1 so that 
deceleration under maximum reverse drive reduces the 
velocity to zero at precisely the command value of the 
output. The drive is then set to zero. 



The ideal relay permits only two conditions; full 
acceleration and full deceleration TRef. 5]. 

From the Bang-Bang control law, we can derive the 
switching function which makes the arror go to zero by using 
the proper switching time. From Newton's second law; 



F = m CRE 

., d . F 

CRE (CRE) 

dt m 



U 



CRE 

But at t 
CRE 

CRE 



CRE dt = 0 t + k1 

0 , ere = 0 and k1 = 0, Therefore 
d 

(CRE) = U t 

dt 

CRE dt = 0 t^ + k2 



(4. 1) 
(4.2) 
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From the equation 4.1 
t = CSE / U 

t^ = (CRE / 0)^ (4.3) 



Substitute equation 4.3 into equation 4.2 

CRE = (0/2) (CRE / 0)^ + k2 

= (CRE)P / (2 0) + k2 (4.4) 

where k2 is integration constant. But if we apply a full 
deceleration at the halfway point, rhe equation 4.4 becomes 

F = (CEE ICREI) / (2 0) + CRE (4.5) 



and is called the ERROR FUNCTION. U will be 
0 = ±G 
or 

0 = -(G) SIGN(F) (4.6) 



Equations 4.5 and 4.6 represem: the SWITCHING FUNCTION which 
makes the error go to zero in the minimum time. The 
switching function and the block diagram of the on-off 
controller are depicted on the Figures 4.1 and 4.2.. And we 
can obtain the cross range error, CRE, by doubly integrating 
D with the initial condition of CRE. We have 



CRE = U dt + CRE(O) 

CRE = CRE dt + CRB(O) (4.7) 

The simulation results of these equations are given on 
Figures 4.3 through 4.7 and the computer program is attached 
(see Appendix B) . 
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Figure 4, 1 Parabolic Switching Function 




Figure 4.2 Block Diagram of ON-OFF Controller 
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Figure 4,3 CRE versus Time 




Figure 4.4 CRE versus Time 
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Figure 4,5 P versus Time 
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Figure 4.6 CHE versus CBE 
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V- MSIC COM MAND TO LI NB- OF- SIGHT SIMOLATION 
A. SCENARIO 

The engagement was designed with the ground tracker and 
missile launching unit located at taa origin. 

The target was flown accross rhe first quadrant from a 
position 4000 meters on the x-axis and 1000 meters on rhe 
y-axis (4000,1000), The velocity vector of the target was 
parallel to the x-axis and magnitude was 250 meters per 
second. 

Since most missiles need a few seconds of boost, the 
missile is not controlled during this time. We assumed that 
the missile was controlled after one second from the firing 
time and controlled by PROGRAMMED 3UIDANCE up to this time. 
After the time of missile "capture'*, the missile was 
controlled by the on-cff, TVC method with the LOS guidance 
law. The simplified flow chart is shown in the Figure 5.1. 

In order to simplify the problem, we assumed that: 

1) the velocity vector of missile , Vm, was parallel to 
the LOS between the target and origin and the magni- 
tude of 7m was constant, 503 meters per second; 

2) the LATAX was applied to the missile at right angles 
to the LOS. This was a reasonable assumption for 
this kind of missile. So the angle <p * in the 
Figure 2.2 is almost same to angle 

3) the measurement noise was zero so we omitted the 
filter , F (s) ; 

4) the magnitude of LATAX was 150 meters /second which 

was about 15 Gs. ' 

The geometry depicted in Figure 5.2 summarizes the geometric 
situation. 
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Figure 5.1 Siaplified Flow Chart of Basic LOS Conaand 
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Figure 5.2 Geoaetry of Basic LOS Guidance 

For mathmatical convenience of simulation, we need to 
define the sign of the CRE and the LATAX as follow; 

♦|CHE| : When the missile position is upper-side of LOS 

-|CRE| : when the missile position is lower-side of LOS 

♦ I U| : when the LATAX is upward direction 

- |0( : when the LATAX is downward direction 

This sign was based on the positive (^ 5 ,^ which is defined 

when is greater than 6 ^ • 

B. PBOGBAHBEO GOI DANCE PHASE 

Since the major emphasis of this paper was on-off 
control, we assumed that the missile flew along the LOS 
during the programmed guidance phase. But, in a practical 
situation, there is some cross range error which is occured 
by disturbances such as wind, propulsion system and 
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autopilot time delay, ere. Hence we made initializarion 
errors, and the on-off control started with these errors. 

C. ON-OFF, THROSr VECTOR, MISSILE CONTROL 

The detail of the on-off control was discussed before, 
hence we applied this to the LOS guidance scheme. The block 
diagram of this system is depicted in the Figure 5.3 
CRef. 3]. 

In order to determine the CRE, the tracker estimates the 
missile's range (Rm) , by the elapsed time of flight and the 
missile's velocity profile. The program of this simulation 
is attached in appendix C. 




Figure 5.3 Block Diagra* of the Basic LOS Command 



D. SiaOLATION RESDLTS 



Figure 5.4 shows the missile and target geometry in X-Y 
plane. The missile intercepted the target at the point 
a(2605,1000) with the almost zero miss distance. 
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Figure 5.5(a) shows the distance between target and 
missile versus time. The distance decreased linearly and 
neared zero at the time at 5.53 seconds. 



— 1 



MSL - TGT TRRJECTORY 
s T-CQOROINRTE VS. X-COGRDINRTE 




Figure 5.4 The Basic LOS Command 

Figure 5.5(b) shows the CRS versus time. The initial 
cross range error at the missile capture time one second 
after firing was about 50 meters. Since the CRE at the 
"missile capture time” was positive, the CRE increased 
initially. So the maximum CRE was about 58. 2 meters at the 
time 1. 330 seconds. Then it decreased to almost zero meter 
at 2.55 seconds. To gat a faster response, we should 
increase the magnitude of the LATilX. We should note here 
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that the CRE does not maintain zero value because we dii nor 
consider the target motion terms in this phase of the simu- 
lation [Bef. 2]. So the missile had some small cross range 
error and the BANG-BANG controller had tried to reduce this 
error in a chatter-mode. 

Figure 5.6(a) shows the CRE versus time. Figure 5.6(b) 
shows the CRE versus CRE. As we expected this curve followed 
the SWITCHING FUNCTION as shown ia figure 4.4. Figures 
5.7(a) and 5.7(b) show the F versus time and the 0 versus 
time . 

This program was tested using maneuvering targets and 
the results were almost same except the impact position. The 
results of this simulation were shown on the Figures 5.8 
through 5.11 and the program is attached in Appendix D. The 
comparison of these simulations is summarized in Table I. 
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Figure 5.5 The Basic LOS Co«iaand 
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Figure 5.6 The Basic LOS Comnand 
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Figure 5.7 The Basic LOS CoBoand 
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Figure 5.8 The Haneuvering Target 




L 



35 



Figure 5.9 The Haneuvering Target 
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Figure 5.10 The Maneuvering Target 




Figure 5.11 0 versus Time for a Haneuvering Target 



TABLE I 

The Basic LOS Command Simulation Eesult 



NON- MANEUV RIHG TARGET 
, timeicontrol^ =1.0 sec 



CRE (0 



9 10(m) 



CRED(O) = 49.832(m/sec) 

. time (max. CR E) 1.33 sec 
■ = 5B. 1 



CHE (max) 



84 (m) 



. time (intercept) = 5.58 sec 
(Xm,Ym) = (2604.7,999. 88) 

CRE = 4.3344 S-5 

distance= 0.34894 (m) 



MANEUVERING TARGET 

1.0 sec 
49.907 (m) 

49.891 (m/sec) 

1 . 33 (sec) 

58.201 (m) 

5.61 3 ec 
(2597. 9, 1057.7) 
-6.5448 S-6 
0.47889 (m) 
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71. PSE DDO- LOS COMHA HD SIMOLATION 



The guidance scheme ci the lead angle command is almost 
the same as that of the basic LOS command. Instead of the 
tracker- to-tar get line-of- sight this guidance scheme uses 
the tracker-to-estimated impact point and is called 
’•synthetic line-of- sight” (SLOS) , or "pseudo line-of-sight". 
The missile is controlled to fly along this pseudo line-of- 
sight. The block diagram of this system is easily modified 
from that of the basic LOS and is shown in Figure 6. 1. 

The estimated impact point at the instantious time is 
calculated by using the "time to go" (Tg) and the "closing 
velocity" (Vc) between the target and the missile. 




Figure 6.1 Block Diagrae of the Pseudo LOS Command System 
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The "closing velocity" and "time to go" are calculated as 
follow; 



simulation result is almost same as in the basic LOS case. 
On other hand, this guidance scheme is poor in a ECM situ- 
ation. In order to compare the results we used the same 
data as that of the basic LOS command. These are shown in 
Figures 6.2 through 6.5 and the summarized results are shown 
in Table II. The computer program is attached in Appendix E. 




Tg = (distance between target and missile) / 7c 




The missile goes to the impacr point directly. The 



TABLE II 

The Pseudo LOS Command Simulation Result 



W iU W I Vw 

CRE(O) 
CRED (0) 



time (control) 



1.0 sec 
49.912 (m) 

49 . 828 (m/sec) 



time (M AZ.CRE) 
CRE (max) 



1.33 sec 
58. 184 (m) 



time (intercept) 
(Xm, Ym) 



5.58 sec 
(2604.7, 999. 

) 



999.92) 

). 0 ) 





miss-distance 



j 
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Figure 6.2 The Pseudo-LOS Command 
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Figure 6.3 The Pseudo-LOS ComiBand 
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Figure 6*4 The Pseudo-LOS Coanand 
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VII. SiaOLA TIOSS WITH TWO^LBVEL RELAY AND SAT ORA TION COraROL 

The LOS guidance with an "ideal" relay has been 
discussed. In this chapter, the effect of the different 
types of lateral acceleration demand are discussed. In 
order to compare the results with the previous simulations, 
the same parameters were used. 

A. TWO-LEVEL RELAY 

The large magnitude of the LATAX makes a fast response. 
But in rhe case of small CRE, a smaller magnitude of LATAX 
is needed. This idea was developed in a "two- level" relay as 
shown on the Figures 7.1(a;b) . The shaded area on Figure 
7.2(b) shows the region of a lower level of LATAX in the 
"CRE verse CRE" phase plain. It provided the minimum over- 
correction. The computer programs were easily modified by 
adding one statement, 

IF ((ICRE I ♦ ICREDI) . LT. M) G = 3/(N1/N2). 

We used the values 150 m/sec^ for N1 and 15 m/sec^ for 
N2 and 1.0 for H in the simulations of the basic LOS command 
and the pseudo LOS command. The results were almost the same 
as the previous, except in the figure for "U versus time". 
Table III summarized these simulation results and Figures 
7.2, 7.3 and 7.4 show the "0 versus time" of each case. The 
programs were attached in Appendix F, G and H. 
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Figure 7.1 Two-Level Belay 



I __________ _________ 

TABLE III 

Two-Level Relay Control Result 



UNIT 


NON-MVR 


MVR-TGT 


PSEUDO-LOS 


. time (control) 
CRE (0) 

CRED (0) 


1.0 

49. 910 
49. 832 


1 . 0 

49.907 

49.891 


1.0 

49.910 

49.828 


.time (MAX. CRE) 
CRE (max) 


1. 33 
58. 184 


1.33 

58.201 


1.33 
58. 184 


. time (intercept) 
Xm 
Ym 
Xt 
Yt 
CRE 

miss-distance 


5.58 
2604.7 
99 9.88 

2605.0 

1000.0 
-2, 63E-6 
0. 34894 


5.61 

2597.9 

1057.7 

2597.5 

1 057. 5 

-3.63E-6 

0.47889 


5.58 

2604.7 

999.92 

2605.0 

1000.0 
5.03E-8 
0. 35137 
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Figure 7.2 Q versus lime for the Basic LOS Guidance 

B. SATORATING LINEAR CONTROL 

In the previous section the two-level relay was 
discussed. The "sautrating linear control" as depicted in 
the Figure 7.5(a) and 7.5(b) was also studied. The shaded 
area on the Figure 7.5(b) shows the region of linear control 
in the "CRE versus CRE" phase plane. The computer programs 
were easily modified by adding one statement, 

IF (ABS (F) . LE. M) U = -G * F / M. 

The value of "N" determines the linear region for F. 
The Figures 7.6 through 7.8 show the simulation results of 
the basic LOS command against the non-maneuvering target 
case for "M" equal 1, 5 and 10. When choosing the value "M" 
equal to "one", the intercept time and miss distance are 
almost the same as the counterpart of the ideal relay case. 
Hence the saturating linear control can be used in practice 



h6 




Figure 7.3 tl verse Tiae against the MVR Target with 

Two-Level Relay 




Figure 7.4 0 verse Tiae of the Pseudo-LOS with 

Two-Level Relay 
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instead of the ideal relay by choosing a proper value of 
"M”. The sumniarized results are in the Table IV. 

Figures 7.9 through 7.11 show the results of the maneuvering 
target case and Figures 7.12 through 7.14 show the results 
of the pseudo-LOS case. These programs are given in 
Appendix I, J and K. 




Figure 7.5 Linear Switching Relay 
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I ABLE 17- 

Saturating Linear Control Result (H=1) 



1 



UNIT 


NON-MVR 


M7R-TGT 


PSEUDO-LOS 


. time (cent ro 1) 
CRE (0) 

CRED (0) 


1.0 

49. 910 
49. 832 


1 . 0 

49.907 

49.891 


1.0 

49.912 
49. 828 


.time (MAX. CRE) 
CRE (max) 


1.33 
58. 184 


1 . 33 

58.201 


1.33 
53. 184 


. time (intercept) 
Xm 
Ym 
Xt 
Yt 
CRE 

miss-distance 


5.58 

2604.7 

999.84 

2605.0 

1000.0 
-.0349 
0. 3507 


5.61 

2598.0 

1057.6 

2597.5 

1057. 5 
-.0702 
0 . 484 1 


5.58 

2604.7 

999.88 

2605.0 

1000.0 
-.0349 
0.3492 



L 



M3L - TGT TRAJECTORY MSL-TGT DISTANCE VS. TIME 

-COORDINATE VS. X-COORDINRTE 3 SATURATION CONTROL 
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Figure 7,6 The Basic LOS Command (Saturation) 
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Figure 7,7 The Basic LOS Command (Saturation) 
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Figure 7-8 The Basic LOS Command (Saturation) 
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Figure 7.9 Haneuvering Target (Saturation) 
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Figure 7.10 Haneuvering Target (Saturation) 
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Figure 7.11 Maneuvering Target (Saturation) 



MSL - TGT TRRJECTORY MSL-TGT DISTANCE V5. T 

-CQOnOINflTE VS. X-COORDlNflTE 3 PSEUDO - LOS COMMAND 





0) 

a 

•H 

tn 

p 

tn 

u 

OJ 

> 

Q) 

0 
C 
(T3 

cn 

•H 

Q 

4J 

(U 

C7> 

}-l 

03 

1 

a 

rH 

•H 

cn 

01 

-H 

s: 






>1 

u 

0 

4J 

0 

CJ 

•f-1 

03 

M 

E-t 

iJ 

01 
CTi 
W 
03 

I 

OJ 

rH 

•H 

W 

03 

•H 



03 



I 

I 

I 

I 



I 

I 



I 



I 

I 



I 



I 



I 



I 



I 



I 

I 



I 



I 

I 



56 



Figure 7.12 Pseado-LOS Command (Saturation) 
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Figure 7.14 Pseudo-LOS Command (Saturation) 



VIII. COHCLDSION 



The comparision of results for the ideal relay, rwo- 
level relay and saturating linear conrrol against the basic 
LOS command and pseudo-LOS command against a non-maneuvering 
target and maneuvering target are provided in the Tables V, 
VI and VII. These simulation results clearly demonstrate 
that "on-off” control of a missile is highly desireable and 
that "saturating linear control" of a missile has little 
adverse effects compared to an "ideal relay" control. 



TABLE V 



Comparison of the Basic LOS Command Results 



UNIT 



IDEAL 

RELAY 



TWO- LEVEL 
RELAY 



SATURATION 

CONTROL 



time (ccntrol) 
CRE (0) 

CRED(O) 



1.0 

49. 910 

49. 832 



1 . 0 

49.910 

49.832 



1.0 

49.910 

49.832 



time (MAX. CRE) 
CRE (max) 




1.33 
58.1 84 




Xm 

Ym 

CRE 



time (intercept) 



5.58 



2604.7 
999. 88 
4.33E-5 
0.34894 



2604.7 
999.88 
4.33E-5 
0 . 34894 



5.58 

2604.7 

999.88 



miss-distance 




.j 
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TABLE VI 








Conparisoa of 


the Maneuvering Target 


Results 




UNIT 


IDEAL 


TWO- LEVEL 


SATURATION 




• •• — 


RELAY 


RELAY 


CONTROL 




. time (ccntro 1) 


1.0 


1 . 0 


1.0 




CRE (0) 


U9. 907 


49.907 


49.907 




CRED(O) 


49. 891 


49.891 


49. 891 




. time (MAX. CRE) 


1.33 


1 . 33 


1.33 




CRE (max) 


58. 201 


58.201 


58. 201 




. time (interc sp t) 


5.6 1 


5.61 


5.61 




Xm 


2597.5 


2597.5 


2597.5 




Ym 


1057.5 


1057. 5 


1057.5 




CRE 


-6. 54E-6 


3 .63E-6 


-0. 0702 




miss-distance 


0.47889 


0 . 47889 


0.4841 


1 








J 








TABLE VII 








ConpaclsoQ of 


the Pseudo-LOS Command 


Results 




ONIT 


IDEAL 


TWO- LEVEL 


SATURATION 




— ” — 


HEL AY 


RELAY 


CONTROL 




. time (contro 1) 


1.0 


1 . 0 


1.0 




CRE (0) 


49. 912 


49.910 


49.912 




CRED(O) 


49. 828 


49.828 


49. 828 




.time (MAX. CRE) 


1.33 


1.33 


1.33 




CRE (max) 


58. 184 


58.1 84 


58. 184 




.time (intercept) 


5.58 


5.58 


5.58 




Xm 


2604.7 


2604.7 


2604.7 




Ym 


999.92 


999.92 


999.92 1 




CRE 


2.96E-5 


5.03E-8 


2. 96E-5 




miss-distance 


0.35137 


0.35137 


0.35137 


1 








j 
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iPPMDIX i 

ViHIABLES LIST 



DIAGRAM 

VARIABLES 


COMPUTER 

VARIABLES 


NOUN DESCRIPTION 


CRE 


CRE 


cross-range- error 


CRE 


CRED 


rate of cross-range-error change 


F 


F 


error function 


G 


G 


magnitude of lateral acceleration 


0 


0 


missile’s lateral acceleration 


dt 


SIGT 


angle between the LOS to target 
and X-axis 




SIGM 


angle between the beam to missile 
and X-axis 


dmt 


SIGMT 


angle difference between m and 


Vc 


VC 


crossing velocity 


tg 


TG 


time to go 


t (control) 


TCON 


beginning time of on-off control 


SLOS 


SLOS 


synthetic line-of-sight 


Xm 


XM 


X-coordiaate of missile position 


Ym 


YM 


Y-coordinate of missile position 


Xt 


XT 


X-coordinate of target position 


Yt 


YT 


Y-coordinate of target position 


Vm 


VM 


velocity of the missile 


Vt 


VT 


velocity of the target 


Vmx 


VMX 


X-component of missile’s velocity 


Vmy 


VMY 


Y-component of missile's velocity 


Vtx 


VTX 


X-component of target's velocity 


Vty 


VTY 


Y-component of target's velocity 
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AP PENDI X B 

PaOGRAH OF THE SHICHING FONCTION 



TITLE 

TITLE 

TITLE 

INTGER 

CONST 

INITIAL 



DERIVATIVE 

NOSORT 

♦ 



BANG-BANG CONTROL 
SWITCHING FUNCTION 
♦ YEDN, J.Y. * 
NPLOr 
NPLOr=1 

CRE = 1.0 
CRED = 0. 

ACRED = 0. 



SAMPLE 



TERMINAL 



G = 1 . 0 
F = CRE 
U = -G * 

CRED 

CRE = INTGRL^CRE, CRED) 
ACRED = ABS(CRED) 

IF(CHE.LE. 0.0) CALL BNDJOB 



+ (CRED*ACRED) /(2*G) 
SIGN(1 . ,F) 

= INTGRL(0.,U) 



CALL 

CALL 

CALL 

CALL 

CALL 



DRWG 

DRWG 

DRHG 

DRWG 

DRWG 



1 , 1, TIME, CRE) 

2.1, CRE,CRED) 

3. 1, TIME,F) 

4.1, TIME, CRED) 

5. 1, TIME,0) 



CONTRL 

PRINT 

END 

STOP 



CALL ENDRW (NPLOT) 

FINTIM=2. 1,DELT=0.01 ,DELS=0.01 
0. 1 ,G, F, 0, CRE, CRED 
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PB06BAS OP 



APPENDIX C 

THE BASIC COttMAMD TO LOS 



TITLE BASIS COMMAND TO L.O.S 

TITLE WITH IDEAL RELAY 

TITLE YEUN, J.Y. 

INTEG RKSFX 

INTGER NPLOr^KILL 

CONST NPLOT=l ,TCON=1.0 

CONST VM=50 0. ,VT=250,PI=3. 1 41 593, KILL=0 

INITIAL 

XTO = 4000. 

YTO = 1000. 

GAMT = PI 
F =0. 

♦ 

DERIVATIVE 

: 4 c:(c 4 t 4 c:tc: 4 c :4(4e:4c4c 4 c 4c4c 4c 4c 4c4c4c4c4>4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c4c 

NOSORT 

4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4t4t4:4c4c 

TARGET PARAMETERS ♦♦♦ 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c4c4c4c4c4c4c4c4c4c 



VTX = VT*COS (GAMT) 

VTY = VT*SINjGAMTi 
XT = VTX*TIME + XTO 
YT = VTY*TIME + YTO 
SIGT = ATAN2(YT,XT) 

IF (TIME.GE.TCON) GO TO 50 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4c 4c4c 4c 4c 

MISSLE PARAMETERS *♦* 

4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4« 4c 4c 4c4c 4c 4c 4c 4«4c 4c 4c 

* PROGRAMMED GUIDANCE ♦ 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4C 4c 4c 4c* 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4= 4c 4c4«4c4c 4c 

♦ 

SIGM = SIGT+0. 1 
VXM = VM * COS (SIGM) 

VYM = VM ♦ SIN (SI3 M) 

XM = INTGRL(0. ,VXM) 

YM = INTGRL(0. ,VYM) 

RM = SQRT(XM**2 + YM**2) 

SIGHT = SIGH - SIGT 
CRE = RM * SIN (SIGHT) 

CRED = DERIV(0. ,CRE) 

GO TO 200 

4e4c# 4: 4e 4c* ******************** ** 

* ON-OFF GUIDANCE (BANG-BANG CONTROL) ********************* 

**♦* *********** *** ********* ** 

* 

50 CONTINUE 

* 



G = 150. 

ACRED = ABSjCRED) 

F = CRE + (CRED*ACRED)/(2*3) 
U = -G ♦ SIGN(1..fL 
CRED = INTGRLIO.O, U) 

CRE = INTGRL(CRE,CRED) 

RM = VM*TIME 
A3 = CRE/RM 
SIGHT = ARSIN(A3) 

SIGM = SIGT + SIGHT 
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XM = aa * cos (sioa) 

YM = RM * SIN(SISa) 

♦ * ♦#♦♦♦♦♦♦* * ♦♦♦ ♦ # 4c4e:tc4e:^4c:te *^5* :(c:4c:4c ^ ^ ^ ;jt ♦* 

200 CONTINUE 

:^:4e34c:4e :4e:4e:4e:4c:4e34c:4e:^:^:4e:$c 3 «c:4c3«c 3^ :«c:^:4(:^34e:4e:4c:4e :#s:#6:C64t:^e:#6«:^::^t3ie:^s:}c3^C3«c:Jc:^c:^s:^:ic:Cs*:«c:#e:Cc:^s«^3jc4c;4c^« 

*** MISSION RESULT **♦ KILL = 0 ; TGT HISSED 

KILL = 1 : TGT DESTROYED 

XDIST = XT-XM 
YDIST = YT-YM 

DIST = SQRT (XDIST**2 + YDIST**2) 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

IF (XM .GT. (XT+30))CALL ENDJOB 

SORT 

OUTPUT AND PLOT CONTROL CARD ♦*♦*♦***♦**♦*#♦♦♦**** 

♦ :«c :(c 4c 4c :)e :<c:4( 4= ♦♦^♦:«'*>S'**^****#1'** ♦♦*=>! ♦♦♦♦♦♦♦**♦ 



SAMPLE 



TERMINAL 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



DRWG 

DRWG 

DRWG 

DRWG 

DRWG 

DRWG 

DRWG 

DRWG 



1 , 1,Xa,YM) 
1,2,XT,YT) 

2. 1, TIME, DIST) 

3 . 1, TIME ,CRE) 

4. 1, TIME,CRED) 

5 .1, CRE,CRED) 



7 , 1,TIME 



6, 1,TIME, F| 
r N) 



CONTRL 

PRINT 

END 

STOP 



CALL ENDRW(NPLOT) 

FINTIM=6. 0, DELT=0.00 1 ,DELS=0. 0 03 
0. 005, XM,YH,XT,YT,CRE,CRSD,DIST,KILL 
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APPM2I1 2 

PHOGRAB OF TBE BANEOVERING TARGET 



TITLE 


BASIC COMBAND TO 


LOS 


TITLE 


WITH BANEUVERING 


IGT 


TITLE 


*** YEUN, J.Y. 


# *# 


INTEG 


RKSFX 




INTGER 


NPLOT,KILL 




CONST 


NPLOT=1 ,TCON=1 . 0 




CONST 


VB=500. ,7T=250, PI 


= 3. 


INITIAL 


XTO = aooo. 
YTO = 1000. 






GABT = PI 
F = 0. 





*♦♦♦*♦♦*#**♦*♦♦ *♦♦*******♦♦ :4c :«(:«( 4c^ :4c :^:4c 4c# 4c 

DERIVATIVE 

4c4c4c*4c4i:tc«**««:^*4c*4c4c4c«4:««4c4c4c:«t 4c4c4t***’)e*******#*#*****4c***** ♦**:««* 

NOSORT 

4c4c*4c »♦♦***♦*** #****♦♦*♦*# 

*** TARGET PARAaETERS ♦** 

4c 4c 4c # 4c * # * 4c « * * 4c « 4c # 4c « 4c * # 4c « ★ « 

VTX = VT*COS <GAMT) 

VTY = VT+SINjGABT) 

XT = VTX*TIME + XrO 

YT = 1 00*SI N (0. 5*PI*TiaE) <-YT0 

SIGT = ATAN2(YT,XT) 

IF (TIME.GE.TCON) GO TO 50 

4; 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4c4c4c 4c 4c4c4c 4«4c 

MISSLE PARAMETERS 

ic :«c :«c :4c :4c :4e 4c « :4c :4c :4c « « >9t :4t :4t « 4e :4c « ^ * :4c :«c 

♦ PROGRAMMED SOIDANCE 

♦ 

SIGH = SIGT+0. 1 
VXM = Va * COS (SIGH 
VYM = va ♦ SIN (SI3 a 
xa = INTGRLfO. ,vxa 
YB = INTGRL(0. .VYa, 

Ra = SQRTfXM**2 ♦ YB**2) 

SIGBT = SIGN - SIGT 
CRE = RB * SIN (SIGBT) 

CRED = DERIV (0. ,CRE) 

GO TO 200 



4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c # 4c 4c 4c4c 4c 4c 4c 4c# 4c 4c 4> 4c 4c4c 

♦ ON-OFF GUIDANCE (BANG-BANG CONTROL) 

4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c 4c 4c 
4c 

50 CONTINUE 

4c 

G = 150. 

ACRED = ABS(CRED) 

F = CRE ♦ (CRED*ACRED)/ (2*G) 
D = -G ♦ SIGN(1..F) 

CRED = INTGRL(0.0, U) 

CRE = INTGRL (CRE, CRED) 

ACRE = ABS (CRE) 



4c 4c 4c 4c 4c 4c 4c 4c# 4c* 4c #4( ## ### *# 



RB •= Va*TIME 
A3 = CRE/RB 
SIGBT = ARSIN(A3) 
SIGN = SIGT + SIG 



BT 
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XM = RM ♦ COS (SIGM) 

YM = RM ♦ SIN (SIGM) 

200 CONTINOS 

♦ 

**:«c:ic:<c3jc*3jc:jc4c:4c:#c:jc:0c* #:{c:0c :«c * * *:«c* :«c:{c :«c :jc :«c * * * :jc :«c :«c :«c :^c * * * :«c :}c * :«c:^:4c4c 3}c4c^:«c* :jc :Jc :jc :4: * 

♦** MISSION RESULT *** KILL = 0 ; TGT MISSED 

:<e:<c RILL = 1 ; TGT DESTROYED 

XDIST = XT-XM 
YDIST = YT-YM 

DIST = SQRT (XDIST**2 + YDIST**2) 

IF (DIST ,LE,5) KILL = 1 

IF (DIST .GT.5[ KILL = 0 

IF (XM .GT. (XT + 30))CALL ENDJOB 

* 

SORT 

4c{c4c:Cc 4c 4c4c4c^:4c4c:4c:tc4c *♦ ♦* 4«***4i4t4t** ♦4t4c * cCt#* ♦ ♦♦ # 4c * ♦ ♦ ♦ #4= # ♦ 4c4c4t 4c4c 4c 4tc^ 4c4c 

4e4c4t4c4e4t4c4c OUTPUT AND PLOT CONTROL CARD 4c4c45 4c4e4c4c4c4c4c4c4c4c4c4c4c4c4c4e4c4c 
4c4c4c4c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4c4c4c 4 c 4 c 4 c 4 c 4 e 4 c 4c4c4c4c4c4c 4c 4c4c4c 4c4c 4c4c4c4c 4e 4c4e 4c 4t 4t 4c 4c4c 4c 4= 4c 4c 4e 4c 4c 4c 4c 4c 4c 4c 4c4c 

SAMP L E 

CALL DRWG(1,1,XM,YM) 

CALL DRNG (1 ,2, XT, IT) 

CALL DEWG (2 , 1,TIME, DIST) 

CALL DRWG (3 , 1,TIME ,CRE) 

CALL DRWG (4 , 1,TIMB,CRED) 

CALL DRWG (5 ,1,CRE, CRED) 

CALL DRWG (6 , 1,TIME, F) 

CALL DRWG (7, 1, TIME ,U) 

TERMINAL 

CALL ENDRW(NPLOT) 

CONTRL FINTIM=5. 9, DELT=0.00 1, DELS=0. 0 03 
PRINT 0.005, XM, Y M, XT, YT,C RE, CRED, DIST, KILL 

END 
STOP 
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APPENDIX B 

PaOGRAH OF THE COMMAND TO PSEUDO-LOS 



TITLE PSEUDO - LOS COMMAND 

TITLE WITH IDEAL RELA? 

TITLE *** YEUN, J.Y. ** * 

INTEG RKSFX 

INTGEE NPLOT,KILL 

CONST NPLOT=1 ,rCON=1.0 

CONST VM=500. , V T=250 , PI= 3 . 1 4 1 5 93 , KILL=0 

INITIAL 

XTO = 4000. 

YTO = 1000. 

TG = 0. 

f = 0. 

GAMT = PI 

♦ 

4 ::(c:(c:tc :(( :(c <c:(c:(c:(c4^:4c:(c:4c *#:jc:(t*#:(c**#*##**#***sOc*^ :(c:(c### 4 c 

DERIVATIVE 

♦ « * :(c:(cic:4c^ ♦**♦*♦*♦♦*♦♦ * ♦4c:{c s(c 4: * ** :4c:(e:(c:<c * 

NOSORT 

*** TARGET PARAMETERS *** 

1 «C« :«c «« :4c # :)c «:<( 4c «:<(««««« 4c :«c :«c :((4c « « 

VTX = VT*COS(GAMT) 

VTY = VT*SIN]GAMT) 

XT = VTX^TIME + XTO 

YT = VTY*TIME + YTO 

SIGT = ATAN2(YT,XT) 

XLOS = XT + ^X*TG 
YLOS = YT + VTY*T3 
SLOS = ATAN2 (YLOS, XLOS) 

IP (TIME.GE.TCON) GO TO 50 

4c4e4c*=)t4c4c4f#4c**4c4c4c*:«t4c4c4«4t4t*4c4c 

MISSLS PARAMETERS 

:4e:4e:4c:4e 4c 3jt 

♦ PROGRAMMED GUIDANCE 4c4c4e4c4c 4c4c4c4c*4c4c4c4c4c4e4c4c4c4c4e4c4c*4c4c4c 4 c4c4c4k4c4c**4c4c* 

♦ 

SIGM = SLOS + 0.1 



VH * COS (SIGM 
VM ♦ SIN(SIGM 
INTGRL(0. ,VM X 
INTGRL (0. eVM Yj 
SQRT(XM**2 ■ 



VMX 
VMY 
XM 

RM = SQRT(XM**2 YM**2) 
SIGMS = SIGM - SLOS 
CRE = RM * SIN (SIGMS) 

CRED = DERIV (0. ,CRE) 

GO TO 2 00 

4c4c4c4c 4c 4c4c4c4c4e4c :<t4c4t 4c 4c 4c4c 4c 4c 4c 4< 4c 4c 4c 4c 4c 4c 4c 

♦ ON-OFF GUIDANCE (BANG-BANG CONTROL) 

4c4c4c4c4c4c4c4c4c4c4c4t4c4c4c4c4c4ci4«4c4c4c4c4c4c4« 4c 4c 
4c 

50 CONTINUE 

4c 

G = 150. 

ACRED = ABS(CRED) 

F = CRE ♦ (CRED*ACRED)/(2*G) 
D = -G * SIGN(1..F) 

CRED = INTGEL(0.0, U) 

CRE = INTGRL (CRE, CRED) 



4c 4t 4c 4c 4c 4c 4c 4c4c4c4c 4c 4c 4c 4c 4c 4c 4« 4c 4c4c 



RM = VM^TIME 
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A3 = CRE/RM 
SIGMS = ARSIN(A3) 

SIGH = SLOS + SIGMS 
VMX = VM * COS (SLOS) 

VMY = VM * SIN(SLOS‘ 

XM = RM * COS(SIGMi 
YM = RM ♦ SIN(SIGMi 

************************ *^* *** ** 

200 CONTINUE 

♦ 

♦ ♦♦♦ ♦ ♦**♦♦♦ ♦♦♦ * ♦♦♦ * ♦ :Cc:^^:4c :4c4c4c ♦♦♦♦*♦*♦#♦ *♦♦♦ ^:«c:4c4c:tE4c:4c4c 

♦** MISSION RESULT *** KILL = 0 ; TGT MISSED 

KILL = 1 ; TGT DESTROYED 

XDIST = XT-XM 
YDISI = YT-YM 

DIST = SORT (XDIST**2 + YDIST**2) 

VC = SQRI ( fVTX-VMX) **2+ (VTY-VMY) #*2) 

TG = DI ST/VC 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

IF Ua .ST. (XT+30))CALL ENDJOB 

♦ 

SORT 

*♦ ♦♦♦ ******** ****************** **** :ic :«c :«c :Cc4< ^ ** 

OUTPUT AND PLOT CONTROL CARD ********************* 



SAMPLE 



TERMINAL 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



DRWG (1 , 1,XM,YM) 

DRHG (1 ,2,XT,YT) 

DRWG (2, 1, TIME, DIST) 
DRHG (3,1, TIME, CRE)_ 
DRWG (4 , 1,TIME,CRED) 
DRWG (5 , 1 ,CRE, CRED) 
DRWG (6, 1,TIME, F) 
DRWG (7 , 1,TIME , U) 



CONTRL 

PRINT 

END 

STOP 



CALL ENDRW(NPLOT) 

FINTIM=5. 9, DELT=0.00 1 ,DELS=0. 00 3 
0. 1 ,TG, XM,YM , XT, YT,CRE, CRED, DIST, KILL 
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APPENDIX P 

PHOGRAH OF THE BASIC COMMAND TO LDS WITH TWO-LEVEL RELAY 



TITLE 


BASIC COMMAND TO 


L .0. S 


TITLE 


WITH TWO LEVEL RELAY 


TITLE 


♦♦♦♦♦ yeun, J.Y. 




INTEG 


RKSFX 




INTGER 


NPLOT,KILL 




CONST 


NPLOT=1 ,rCO N=1. 0 




CONST 

INITIAL 


VM=500 . ,VT=250,PI 


=3, 141593, KILL=0 




XTO = UOOO. 
YTO = 1000. 
GAMT = PI 





F = 0. 

♦ 

* :4c ♦#:4e**#** ***:4e*^ ♦ :je:4c + aic :4c :)t :4c # 

DERIVATIVE 

♦ ♦ 4c:4c ♦ ♦♦♦*♦* :4c ♦ :4c ♦♦:4c* ♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦ ♦ ♦♦♦♦♦♦♦♦ ♦♦♦ 

NOSORT 

****♦♦♦*♦♦♦♦♦♦*♦*♦**♦♦♦♦♦ 

♦** TARGET PARAMETERS *** 

♦♦♦♦ ♦ ♦♦♦♦♦♦♦♦♦ ♦♦ ♦♦♦ ♦♦♦♦♦♦ 

VTX = VT*COS(GAMT) 

VTY = VT*SINjGAMT) 

XT = VTX*TIME + XTO 
YT = VTY*TIME ♦ YTO 
SIGT = ArAN2(YT-XT) 

IF (TIME. GE .ICON) GO TO 50 
♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦ 

*** MISSLE PARAMETERS *** 

♦♦♦♦ ♦ ♦♦♦♦♦♦♦♦♦ ♦ ♦ ♦♦♦ ♦ ♦♦♦♦ ♦ 

* PROGRAMMED GUIDANCE ♦♦*♦♦**=«'***♦*****♦♦♦*****♦******♦***♦* 

* 

SIGM = SIGT + 0.1 
VXM = VM * COSISIGM) 

VYM = VM ♦ SIN(SI3Mi 
XM = INTGRL(0. ,VXMi 
YM = INTGRL (0. .VY M) 

RM = SQRT(XM**2 + YM**2) 

SIGMT = SIGM - SIGT 
CRE = RM * SIN (SIGMT) 

CRED = DERIV(0. ,CRE) 

GO TO 2 00 

♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦ ♦ ♦♦♦ ♦♦♦♦♦♦♦♦♦ ♦♦ 

♦ ON-OFF GOIDANCE (BANG-BANG CONTROL) 

♦ 

50 CONTINUE 



G = 150. 

ACRE = ABS(CRE) 

ACRED = ABS (CRED) 

IF ((ACRE* ACRED) .LT. 1.) G = 15 
F = CRE + (CRED*ACRED)/ (2*G) 

U = -G * SIGN(1 .,F) 

CRED = INTGEL(0.0, U) 

CRE = INTGRL (CRE, CRED) 

RM = VM*TIME 
A3 = CRE/RM 
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SIGHT = ARSIN(A3) 

SIGH = SIGT + SIGHT 
XH = RH * COS (SIS H) 

YH = RH * SIN (SIGH) 

:<c4c3te:(c ♦ * :{c 4c:{c :{c « « :)c :{c 4c ^ * ^cjc^ic ♦❖♦4: ♦ *♦♦♦♦♦♦ ♦♦ 



200 



CONTINUE 



♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦ ♦ ♦ ♦♦ ♦♦♦♦♦♦ ♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦ ♦♦ 



*** HISSION RESULT *** KILL = 
♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ KILL = 
XDIST = XT-XH 
YDIST = YT-YH 
DIST = SQRT(XDIST**2 



0 ; TGT HISSED 

1 : TGT DESTROYED 



IF (DIST .LE.5 
IF (DIST .GT.5 
IF (XH .GT. (X 



I 



YDIST**2) 
KILL = 1 
KILL = 0 
30) ) CALL ENDJOB 



SORT 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦ ♦♦♦♦♦♦♦♦ ♦♦ 

♦♦♦♦♦♦♦♦ OUTPUT AND PLOT CONTROL CARD *♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
♦♦♦ ♦ ♦♦♦♦♦♦♦ ♦♦♦ ♦ ♦♦♦ ♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦ ♦♦♦♦♦♦♦♦♦ ♦♦♦ ♦♦ 



sample 

CALL DRWG (1 ,1,XH,YH) 

CALL DRWG (1 ,2,XT,YT) 

CALL DRWG (2 , 1,TIHE ,DIST) 

CALL DRWG (3 , 1,TIHE,C3E) 

CALL DRWG (4 ,1,TIHE ,CRED) 

CALL DRWG (5 , 1,CRE, CRED) 

CALL DRWG (6 ,1,TIHE ,F) 

CALL DRWG (7, 1,TIHE,U) 

TERHINAL 

CALL ENDRWiNPLOT) 

CONTRL FINTIH=5. 9, DELT=0.00 1 ,DSLS=0. 0 03 
PRINT 0.005, XH,YH, XT, YT, C RE, CRED, D 1ST , KILL 

END 
STOP 
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aPP^DlX G 

PHOGRAa OF THE aANEOVERING TARGET WITH TWO-LE7EL RELAY 



TITLE 


BASIC COMMAND TO 


LOS 


TITLE 


WITH MANEUVERING 


TGT 


TITLE 


4C4C4C YEUN, J.Y. 


4c 4c4c 


INTEG 


RKSFX 




INTGER 


NPLOT,KILL 




CONST 


NPLOT= 1 ,TCO N=1 . 0 




CONST 

INITIAL 


VM=500. ,VT=250,PI 


= 3. 




XTO = 4000. 
YTO = 1000. 
GAMT = PI 





F = 0. 

:4e:(: :4c :4c 4c 4c 4c 4c 4c 4c 4c 4 * 4c 4c 4c4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c 4c 4c 4c 4c 4c 4c 4c4c 4c4c4c 

DERIVATIVE 

4c4c4c4c 4c4c4c4c4c4c4c4c4c4c4c 4c4c4c4c4c 4c 4c 4c 4c 4< 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

NOSORT 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4c4c 4c 4c 4c 4c 4c 4c 4c 

♦** TARGET PARAMETERS *** 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4c 4i 4c 4c 4c 

VTX = VT*COS <GAt1T) 

VTY = VT*SINiGAMTj 
XT = VTX*TIME + XTO 
YT = 100*SIN (0.5*PI*TIME) +YT0 
SIGT = ArAN2{YT,XT) 

IF (TIME.GE.TCON) GO TO 50 

4c4c4c4c 4c4c4c4c4c4c4c4t4c4c4c4e 4c 4c 4c 4c 4c 4c4c 4 : 4c 

MISSLE PARAMETERS *** 

* PROGRAMMED GUIDANCE 

♦ 

SIGM = SIGT+0.1 

VM » COS (SIGH) 

VM * SIN(SIGMi 
INTGRL(0. ,VXM) 

INTGRL (0. ,VY M) 



VXM 
VYM 
XM 
YM 
RM 
SIGHT 
CRE = 
CRED 



= SQRT{XM**5 *■ YM=«‘*2) 
T = SIGH - SIGT 



RH * SIN (SIGHT) 

= DERIV (0- ,CR E) 

GO TO 2 00 

4c4c4c4c ♦ 4t 2®c4e4t 4 c 4c ♦♦***3Jc* 4c 4c 

* ON-OFF GUIDANCE (BANG-BANG CONTROL) 

4c4c4c4c 4c4c4c4c4c4c4c4c4c4c 4c4c4c4c i4c4(4c4c4c4c4c4c 4c 4c 
4c 

50 CONTINUE 

4c 

G = 150. 

ACRE = ABS(CRE) 

ACRED = ABS (CRED) 
IF((ACRE+ACHED) .LT. 1.) G = 15. 
F = CRE + (CRED*ACRED)/ (2*G) 

0 = -G * SIGN(1..FL 
CRED = INTGRL(0.0, 0) 

CRE = INTGRL(CRE,CRED) 

ACRE = ABS (CRE) 

4c 

RM = VM*TIME 
A3 = CRE/RH 



4c 4c 4c 4c 4c 4c 4c 4c4c4c4c 4c 4c 4c 4c 4c 4c4c4c4c4c 
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SIGHI = iRSIN(A3) 

SIGM = SIGT + SIGHT 
XM = RM * COS (SIGM) 

YM = RM * SIN (SIGM) 

200 CONTINUE 

♦ 

**♦♦♦♦*♦*♦♦ ♦♦♦♦ ♦ :<(:4(4c4c;4c4e:tc4E4c :«(:tc4e 4E4c3tt *4c4c:<c *♦*♦♦♦♦* ♦♦ 



MISSED 

DESTROYED 



*** MISSION RESULT **♦ KILL = 0 ; TGT 

♦♦♦**#♦♦***♦♦♦♦♦♦♦♦*♦♦ KILL = 1 ; TGT 

XDIST = XT-XM 
YDISr = YT-YM 

DIST = SORT (XDIST**2 + YDIST**2) 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

IF (XM .GT. (XT+30))CALL ENDJOB 

^ . 

SORT 

*♦♦♦**♦ *♦♦♦♦**♦♦♦*♦♦♦♦**♦**♦♦* ♦**:«t*:jt** 

OUTPUT AND PLOT CONTROL CARD 

4c:(e:te:4c ^ ^ :fc4c ^ 

SAMPLE 



TERMINAL 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



DRWG 

DRHG 

DRWG 

DHMG 

DRWG 

DRWG 

DRWG 



1 ,1,XM,YM) 

1 ,2,XT,YT) 

2 .1, TIME ,DIST) 

3 . 1, TIME,CRE) 

4 . 1, TIME .CRED) 

5. 1, CRE, CRED) 

6 . 1, TIME ,F) 



DRWG (7 , 1,TIME,U' 



CONTRL 

PRINT 

END 

STOP 



CALL ENDRW(NPLOT) 

FINTIM=5. 9,DELT=0.00 1 ,DELS=0. 003 
0.005, XM,YM, XT,YT, CRE,CRED,DIST, KILL 
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APPENMX H 

PROGBAH OF THE COHHAND TO PSEODO-L3S ilTH TWO-LEVEL RELAY 



TITLE PSEODO - LCS COMMAND 

TITLE WITH TWO-LEVEL RELAY 

TITLE YEON, J.I. ♦♦♦*♦ 

INTEG RKSFX 

INTGER NPLOr,KILL 

CONST NPLOT=1 ,TCON=1 .0 

CONST VM=530. ,VT= 250, PI= 3. 141 593, KILL=0 

INITIAL 

XTO = 4000. 

YTO = 1000. 

TG = 0. 

GAMT = PI 
F = 0. 

♦ 

♦ ♦♦♦ ♦♦♦♦♦♦♦♦♦♦ ♦ ♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦ 45 41^ ^ 

DERIVATIVE 

««*««««««««« *********** 

NOSORT 

*****i)i********** ********* 

*** TARGET PARAMETERS **♦ 

:«e«:4c :$c:4c« :4c « :4c :4c:4c :«( 



VTX = VT*COS(GAMT) 

VTY = VT*SINjGAMrj 
XT = VTX*TIME + XTO 
YT = VTY*TIME ♦ YTO 
SIGT = ATAN2(YT,XT) 

XLOS = XT + VTX*TS 
YLOS = YT + VTY*TG 
SLOS = ATAN2 (YLOS, XLOS) 

IF (TIME. GE. ICON) GO TO 50 
**** ********************* 



HISSLE FARAMErEBS 

4E:4e:4e * 4c 4c :4c4(«4e « « ^:4c:4e :4c :4c 

* PROGRAMMED GUIDANCE 4c4c4c4c4c 4e4c4e4c4c4e4c4c4c4c4c4c4c4c4c4c4e4c4e4c4c4e4c4c4c4c4c4e4c4c4c4c4c 
4c 

SIGH = SLOS + 0.1 
VMX = VM * COS (SIGM) 

VMY = VM * SIN (SIGM I 
XM = INTGRL (0. ,VMXi 
YM = INTGRL(0. .VMY) 

RM = SQRT(XM**2 ♦ YM**2) 

SIGMS = SIGM - SLOS 
CRE = RM * SIN (SIGMS) 

CRED = DERIV(0. ,CRE) 

GO TO 200 

*************************** ** 

* ON-OFF GUIDANCE (BANG-BANG CONTROL) ********************* 

****************** it******** ** 

* 



50 CONTINUE 

G = 150. 

ACRE = ABS(CRE) 

ACRED = ABS(CRED) 

IF ((ACRE+ ACRED) .LT. 1.) G = 15. 
F = CRE + (CRED*ACRED)/ (2*3) 

U = -G * SIGN(1..F) 

CRED = INTGRLIO.O, U) 

CRE = INTGRL (CRE, CRED) 
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I 



RH = VM*TIME 
A3 = CRE/RM 
SIGHS = ARSIN(A3) 

SIGH = SLOS + SIGHS 
VMX = VH ♦ COS (SLOS) 

VHY = 7H * SIN(SLOSi 

XM = RH * cos(sigm;i 

YH = EH * SIN (SIGH I 
200 CONTINUE 

★ 

:Jc:Sc:Jc* :#c ^ jJc :jc :jc ♦ :#c :jc * :jc*^ *:#c :^c :(e 3(e 3$e :$c :$e :Jc * * ;je :«c * ♦ :Jc ajc * :J6 ♦*:0c:Jc :4c :#e:Sc :Jc 

*** HISSION RESULT *** KILL = 0 ; TGT HISSED 

^i*:iti*M**MM*M*M*M*iic****iic KILL = 1 ; TGT DESTROYED 

XDIST = XT-XH 
YDIST = YT-YH 

DIST = SQRT(XDIST**2 + YDIST**2) 

VC = SQRI ( (VTX-VMX) **2+ (VTY-VHY) **2) 

TG = DI ST/VC 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

IF (XH .GT. (XT+30))CALL ENDJOB 



SORT 

*************** ******** ******************************** 
OUTPUT AND PLOT CONTROL CARD ********************* 

Hfiti** ***************************************** ************** 

SAHPLE 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



DR MG 
DR MG 
DBMG 
DR MG 
DRMG 
DR MG 
DRMG 
DRWG 



n,1,XH,YH 
(1 ,2,XT,YT) 

2. 1, TIME, DIST) 

3 . 1, TIME ,CRE) 

4. 1, TIME,CRED) 

5 .1, CRE,CRSD) 

6.1, TIME,F' 

7.1, TIME ■■ 



TERMINAL 

CALL ENDRM(NPLOT) 
CONTRL FINTIM=5. 7,DELT=0.00 1 



,U) 



DELS=0 . 003 



PRINT 

END 

STOP 



0. 1,TG, XH, YH ,XT,YT, CRE ,CRED, DIST , KILL 
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&P PEN DIX I 

OF THE BASIC COMHAHD TO LOS WITH SATORAnOH CONTROL 



BASIS COMMAND TO L.O.S 
MISSLE CONTROL 
WITH SATO RATION CONTROL 
RKSFX 

NPLOr, KILL. CUR 
NPLOT=1 ,ICON=1.0,CDR=1 
VM=50 0 . ,7T=250, PI=3. 141 5 93, KILL=0 
M = 10. 

XTO = 4000. 

YTO = 1000. 

GAMT = PI 
F = 0. 

: 4 c 4 c 4 c: 4 c :(c 4c4c 4c :(c :(c :<(:)[ 4 c 4 c 4 c 4 c 4c4c4c4c4c4c4c4c4c 4c*#:4c4c4c4(3tc4c4e4‘’j‘###4c 4c4c 4c4c4c4< 4c4c4c4(4c* ** 4c4c4e 

DERIVATIVE 

4c 4c 4 c 4 c 4 c 4c4c 4c 4c4c4c 4 c 4 c 4 c 4 c 4c 4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c:<c4c4c4c4c4c:(c4c4c4c4c4c4c4c4c4c4c 4c4c4c4c 4c 4c 4c 4c 4c 4c4c4c 4c4c4c 

NOSORT 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c etc 4c 4c 4c 4c 4c « 

4C4C4C target parameters *♦* 

4 c 4 c 4 c 4 c 4c 4c 4c 4c 4i 4c 4c 4c4c4c 4c 4c 4c4c 4c 4c 4c4c4c4c4c 

VTX = VT*C0S(GAMT) 

VTY = VT*SINjGAMr) 

XT = VTX*TIHE + XTO 
YT = VTY*TIME + YTO 
SIGT = ATAN2(YT,Xr) 

IF (TIME. GE. ICON) GO TO 50 

MISSLE PARAMETERS 

♦ :^e :?c :0c :(c :^c ♦ :0c :0c:0e :0c ♦ :0c :0c :0c :0c 

♦ PROGRAMMED GUIDANCE ♦=0c:Oc*4c ***:0c* + 4c:0c^*+ + :0c:0c4c4c*:0c:0c*:0c* :0c:0c*:0c:0c:0c**:0c:0c:0c 

♦ 

SIGM = SIGT + 0.1 
VXM = VM * COS (SIGM) 

VYM = VM * SIN (SIGM 
XM = INTGRL (0. ,VX M) 

YM = INTGRL(0. .VYM) 

RM = SQRT(XM**2 + YM**2) 

SIGHT = SIGM - SIGT 
CRE = RM * SIN (SI3 MT) 

CRED = DERIV(0. ,CRE) 

GO TO 2 00 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4c 4c 4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

♦ ON-OFF GUIDANCE (BANG-BANG CONTROL) 4c4=4c4c 4c4c4c 4c4c4c4c 4c4c4c4c4c 4c4c4c 4c4c 
4c4c4c4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c4c 4 c 4 c 4 c 4 c 4 c 4 c 4c4c4c 4c 4c 
4c 

50 CONTINUE 

4c 

G = 150. 

ACRED = ABS (CRED) 

F = CRE + f CHED*ACRED) / (2*G) 

U = -G * SIGN(1 .,P) 

IF (ABS(F) .LT. M) 0= -G*F/M 
CRED = INTGRL(0.0, U[ 

CRE = INTGRL(CRE,CRED) 

RM = VM*TIME 
A3 = CRE/RM 
SIGHT = ARSIN(A3) 



PROGRAM 



TITLE 

TITLE 

TITLE 

INTEG 

INTGER 

CONST 

CONST 

PARAM 

INITIAL 
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SIGM = SIGT + SIGHT 
XM = RM * COS (SIGM) 

YM = RM * SIN(SIGM) 

♦ ************* 

200 CONTINUE 

♦ 

********** He**** **He ********* ******************************** 

MISSION RESULT **♦ KILL = 0 ; TGT MISSED 

********************** kill = 1 ; TGT DESTROYED 

XDIST = XT-XM 
YDISr = YT-YM 

DIST = SQRT (XDIST**2 + YDtST**2) 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

SORT 

**** *********************** ******************************** 

******** OUTPUT AND PLOT CONTROL CARD ********************* 
*** * ********** * *** ************ *************** ******* ** *** ** 

SAMPLE 

CALL DRWG (1 ,CUR,XM ,YM) 

CALL DRWG (2 ,CUR,TI ME,DIST) 

CALL DEWG (3 ,CUR,TI ME, CRE) 

CALL DRWG (a ,CUR,CRE,CRED) 

CALL DRMG (5 ,CUR,TIME, F) 

CALL DRMG U »CUR,TIME, U) 

TERMINAL 

IF (CUR .EQ. 3) CALL ENDRW(NPLOT) 

CUR = CUR ♦ 1 

CONTRL FINTIM=5. 65 , D ELT=0. 0 0 1 , DE LS=0 . 00 3 
PRINT 0.005,XM,YM, XT, YT, CRE, CRED, DIST, KILL 
END 

PARAM M = 5. 

END 

PARAM M = 1. 

END 

STOP 
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APPENDIX J 

PR0GBA8 OF THE MANEDVEEING TARGET IITH SATURATION CONTROL 



TITLE BASIC COMMAND TO L 
TITLE (MANEUVERING IG 

TITLE WITH SATURATION CO 

INTEG RKSFX 

INTGER NPLOr, KILL, CUR 

CONST NPLOT=1 ,ICON=1. O.C 

CONST VM=500. ,VT=250,PI= 

PARAM M = 10. 

INITIAL 

XTO = 4000. 

YTO = 1 000. 

GAMT = PI 

F = 0. 



♦ 



OS 

N^ROL 



UR= 1 

3. 141593, KILL=0 



:«c :(c :{c # :«c ♦♦ 



DERIVATIVE 

NOSORT 

:<c:4c«:4c :(e#:0c*««:4e:4c#« « :«c #4c :4c « :<c:4c4c :«c :4c 

TARGET PARAMETERS *** 

*♦♦♦♦♦♦* ♦:4c :4c^* ♦ ♦ 

VTX = VT*COS (GAMT) 

VTY = VT*SIN]GAMTi 
XT = VTX*TIME + XTO 
YT = 100*SIN (0. 5*PI*TIME) +YT0 
SIGT = ATAN2jYT,XT) 

IF (TIME, GE. ICON) GO TO 50 

:4c^4c^ :4c :4c^^#:4( :4c ♦:4c:4c :4c ♦♦♦♦♦♦♦♦♦♦ 

MISSLE PARAMETERS *** 

♦ 3«c ♦ # :«c ♦ ♦ ♦ :^c # ♦ * :|£ ijc 4c ♦ ♦ ♦ ♦ 3^c 

♦ PROGRAMMED GUIDANCE 

♦ 

SIGM = SIGT+0. 1 
VXM = VM * COS (SIGM) 

VYM = VM * SIN (SIGM) 

XM = INTGRL(0. ,VXM) 

YM = INTGRL (0, .VYM) 

RM = SQHT(XM**2 ♦ YM**2) 

SIGMT = SIGM - SIGT 
CRE = RM * SIN (SIGMT) 

CRED = DERIV (0. ,CRE) 

GO TO 2 00 

♦ ♦♦♦♦♦**:4c***»* ♦♦♦♦♦♦ ♦♦♦♦♦♦♦ ♦♦ 

♦ ON-OFF GUIDANCE (BANG-BANG CONTROL) ♦♦♦♦♦♦♦*♦♦♦♦♦*♦♦♦♦♦♦♦ 

♦:4e********^ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦ 

♦ 

50 CONTINUE 

* 

G = 150. 

ACRED = ABS(CRED) 

F = CRE + (CRED^ACRED)/ (2*G) 

D = -G * SIGN(1.,F) 

IF (ABS (F) . LT.M) U=-G*F/M 
CRED = INTGRLfO.O, U)_ 

CRE = INTGRL (CRE,CRED) 

ACRE = ABS (CRE) 

RM = VM*TIME 
A3 = CRE/RM 
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SIGHT = 4RSIN(A3) 

SIGH = SIGT + SIGHT 
XM = RM * COS(SIGM) 

YM = RM » SIN (SIGH) 

^ :(c:«c^:4c4c4e:>; ♦ jjt jjc 

200 CONTINUE 

♦ 

♦ ♦♦♦♦♦♦*«♦♦♦♦♦ ♦♦#******:«c*#* ******* **^^** ****** ********** 

MISSION RESULT KILL = 0 ; TGT MISSED 

*^*****4i*****:i^*:^*****-if kill = 1 ; TGT DESTROYED 

XDIST = XT-XM 
YDIST = YT-YM 

DIST = SQRT (XDIST**2 + YDrST*«2) 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

SORT 

♦ ^i^f:tii:ic*:lic**:lic**:^*ifif:tf ****** ********** 

OUTPUT AND PLOT CONTROL CARD ********************* 

*****:lf:tli*** **** ♦ :4[:tcXe jjc :(c :<c « « :<c ♦♦ ****^*** ** 

SAMP L E 

CALL DRWG (1 ,CUR,XM,YM) 

CALL DRWG (2 ,CUR,TI ME, DIST) 

CALL DRUG (3 , CUR, TIME, CRE) 

CALL DRWG (4 ,CUR,CRE,CRED) 

CALL DRMG (5,CUR,TIME, F) 

CALL DRWG (6 ,CUR,TI ME, U) 

TERMINAL 

IF (COR .BQ.3) CALL ENDRW(NPLOT) 

COR = CUR + 1 

CONTRL FINTIM=5. 65 , D ELT=0. 0 0 1 , DE LS=0 . 00 3 
PRINT 0.005,XM,YM, XT,YT, CRE,CRED, DIST, KILL 
END 

PARAM M = 5. 

END 

PARAM M = 1. 

END 

STOP 
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APPENPII 5 

PBOGRAM OF THE COHMAND TO PSSDDO-LOS WITH SATOHATION CONTROL 



TITLE PSEODO - LOS COMMAND 

TITLE WITH SAIHRATION CONTROL 

TITLE YEUN, J.Y. ***** 

INTEG RKSFX 

INTGER NPLOT, KILL, CUR 

CONST NPLOT=1 ,TCON=1.0,CUR=1 

CONST VM=50 0. ,VT= 250, PI = 3 . 1 41 593 , KILL=0 

PARAM M = 10. 

INITIAL 

XTO = 4000. 

YTO = 1000. 

TG = 0. 

GAMT = PI 
F =0. 

♦ 

DERIVATIVE 

:4c4i:tc :(c 4c4c# ♦ :4c :4c:tc ♦ * :4c****** ****** *** ******* ** **** ******** *** 

NOSORT 

**** ********** ****** ***** 

TARGET PARAMETERS *** 

************************* 



VTX = VT*COS(GAMT) 

VTY = VT*SINjGAMT) 

XT = VTX+TIME + XTO 
YT = VTY+TIME YTO 
SIGT = ATAN2(YT,XT) 

XLOS = XT + VTX*T3 
YLOS = YT ♦ VTY^TG 
SLOS = ATAN2 (YLOS, XLOS) 

IF (TIME. GE.TCON) GO TO 50 
************************* 



*♦* MISSLE PARAMETERS *** 

************************* 

* PROGRAMMED GUIDANCE *******♦****♦**♦*♦♦♦♦**♦*♦♦ ♦♦♦ 

♦ 

SIGM = SLOS + 0.1 
VMX = VM ♦ COS (SIGM) 

— TT M CXXT / CT ^ M , 

XM 
YM 
RM 

SIGHS = SIGM - SLOS 
CRE = RM * SIN(SIGMS) 

CRED = DERIV(0. ,CRE) 

GO TO 2 00 

*** * * ********* * * *********** ** 



VM * SIN (SI3 M 
INTGRL(0, ,VMX 
INTGEL(0. .VM Y 
SQRT(XM**2 ♦ 



YM**2) 



♦ ON-OFF GUIDANCE (BANG-BANG CONTROL) *♦♦♦♦**♦♦♦♦♦♦**♦*♦**♦ 
**** * ****** *** * * ** ********* ** 

♦ 

50 CONTINUE 

* 

G = 150. 

ACRED = ABS (CRED) 

F = CRE ♦ (CRED+ACRED)/ (2*3) 

U = -G * SIGN(1.,F) 

IF (ABS(F)^ .LT.M) U = -G * F /M 
CRED = INTGRL(0.0, U)_ 

CRE = INTGRL (CRE,CRED) 
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RM = 


7 M^TIME 




A3 = 


CRE/RM 




SIGMS 


= ARSIN(A3) 


SIGM 


= SLOS + SIGMS 


VMX = 


V M ♦ COS < 


rSLOS) 


VMY = 


VM * SIN( 


SLOS 


XM = 


R H ♦ COS ( 


SIGM 


YM = 


RM ♦ SIN ( 


SIS M' 



:^:(e :0c :4c :0c :ji:0c * 4c30c ♦ :0c:0e:0c:0c4c:0c:0c:0c:0c:0c:0c:0c:4c:0c:0c:0c:0c:0c:0c:0c:0c:0c ^:0c:0:*:0c :0c*:0s :0c:0c 

200 CONTINUE 

♦ 

:0c :0c :0c 4: :0c :0c :0c :0c :0c :^e :0c ::0c:0c:0c ::0c :0c :0c:0c :0c * * :0c :0c :0c * :0c :0c ** :0c :0c :0c :0c * :0c :0c * :0c :0c :0c :0c :0c * :0c :0c 4c=0c:0c:0c :0c :0c :0c :0c :0c :0c:0c:0c :0c:0e 

*♦* MISSION RESULT ♦♦♦ KILL = 0 ; TGT MISSED 

KILL = 1 ; TGT DESTROYED 

XDISr = XT-XM 
YDIST = YT-YM 

DIST = SQRT (XDIST**2 + YDIST**2) 

VC = SORT ( ( VTX-VMX) **2+ (VTY-VMY) **2) 

TG = DISI/VC 

IF (DIST .LE.5) KILL = 1 

IF (DIST .GT.5) KILL = 0 

SORT 

:t^:t(:4c:4c **♦ ♦*♦♦♦*♦♦ ## ♦ ♦♦♦#♦#*# ♦# 

OUTPUT AND PLOT CONTROL CARD ********************* 

**** 4 c 4 : 4 ; > 4 c:{c:^ 4 c 4 ;:(c 4 ; 4c 4c 4c 4= 4c ♦ 4c 4c 4t 4= 4c 4c 4= 4: 4c 4c 4c 4c 4c 4c 4c 4e 4c 4e 4c 4c 4:4c 4c 4c4c 4c4c4c4c 4 c 4c 4c 4c 4e 4c 4c 4c 4c4c 

SAMP L E 

CALL DRUG (1 , CUR, XM ,YM) 

CALL DRHG (2 ,CUR,TIME,DIST) 

CALL DRUG (3 , CUR, TIME, CRE) 

CALL DRMG (4 ,CDR,CEE,CRED) 

CALL DRWG (5 ,CUR,TIME, F) 

CALL DRWG (6 ,CUR,TI ME, U) 

TERMINAL 

IF (CUR .EQ.3) CALL ENDRW(NPLOT) 

CUR = CUR + 1 

CONTRL FINTIM=5. 65,DELT=0.001 ,DELS=0.003 
PRINT 0.005,TG,XM, YM,XT, YT,CRE, CRED, DIST, KILL 
END 

PARAM M = 5. 

END 

PARAM M = 1. 

END 

STOP 
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